Ultra broadband flat dispersion tailoring on reversed-rib chalcogenide glass waveguide
Zhai Yanfen, Qi Renduo, Yuan Chenzhi, Zhang Wei†, , Huang Yidong
Tsinghua National Laboratory for Information Science and Technology, Department of Electronic Engineering, Tsinghua University, Beijing 100084, China

 

† Corresponding author. E-mail: zwei@tsinghua.edu.cn

Project supported by the National Basic Research Program of China (Grant Nos. 2013CB328700 and 2011CBA00303) and the National Natural Science Foundation of China (Grant Nos. 61575102 and 61321004).

Abstract
Abstract

In this paper, we introduce a horizontal slot in the reversed-rib chalcogenide glass waveguide to tailor its dispersion characteristics. The waveguide exhibits a flat and low dispersion over a wavelength range of 1080 nm, in which the dispersion fluctuates between −10.6 ps·nm−1·km−1 and +11.14 ps·nm−1·km−1. The dispersion tailoring effect is due to the mode field transfer from the reversed-rib waveguide to the slot with the increase of wavelength, which results in the extension of the low dispersion band. Moreover, the nonlinear coefficient and the phase-matching condition of the four-wave mixing process in this waveguide are studied, showing that the waveguide has great potential in nonlinear optical applications over a wide wavelength range.

1. Introduction

Chromatic dispersion is an essential property of nonlinear waveguides, which plays a critical role in many nonlinear processes, such as soliton formation,[16] super-continuum generation,[3,710] four-wave-mixing (FWM) based amplification,[11] and wavelength conversion.[1216] Low and flat dispersion over a broad wavelength range is preferred for optimizing the phase matching condition.[17] Recently, chalcogenide glass (ChG) waveguides have drawn much attention as promising candidates for integrated nonlinear photonic devices.[1829] Many ChG materials have high Kerr nonlinear indexes (n2) and low two-photon absorption coefficients;[3032] hence, they have high figure of merit (FOM) of the third-order nonlinear optical materials.[30,31] However, ChG materials usually have large negative dispersion at the near-infrared and mid-infrared bands, which requires a special waveguide design to compensate for the material dispersion by the waveguide dispersion.[17,22,23,3337]

Strip and rib waveguides are the most popular structures of ChG waveguides. The geometry parameters of these structures can be adjusted to tailor the waveguide dispersions. However, the bandwidth of low dispersion region is always limited. The introduction of a slot structure into the waveguide has been proven as an effective way to tailor the waveguide dispersion.[12,23,3842] Recently, a strip/slot hybrid As2S3 waveguide design has been proposed, which has a flat dispersion band of 249 nm and whose dispersion is limited in ±170 ps·nm−1·km−1.[43] Although this design shows impressive dispersion characteristics as nonlinear waveguides, many difficulties can be expected in its fabrication. On the one hand, the commercial (NH4OH)-based developer using in the UV-lithography or electron beam lithography (EBL) may lead to pinholes and film peeling in ChGs[4345] which impacts the quality of ChG waveguides. On the other hand, both the wet-etching and dry-etching processes on the ChG films are not easy due to the fragility of ChG materials.[27,4549] Recently, we proposed a reverse ridge/slot hybrid As2S3 waveguide structure with an ultra-flat dispersion profile.[23] This structure can be fabricated by the micro-trench filling technique,[50] which avoids the photo lithography process, wet or dry etching processes and lift-off processes on the As2S3 layer. However, to realize this waveguide structure, vertical slot structures of tens of nanometers should be fabricated by dry etching on the silica substrate, and the requirement of the accuracy of the structure parameters is very high to realize dispersion flattening, which is difficult in waveguide fabrication.

In this paper, we introduce a horizontal slot in a reversed-rib As2S3 waveguide. This waveguide can also be fabricated by the micro-trench filling technique; however, the slot in the horizon direction replaces the vertical slots in the previous work[23] to tailor its dispersion property. Hence, the difficulties in the fabrication of the small slot structures in the silica substrate are avoided, leading to a simpler fabrication process. On the other hand, the dispersion of quasi-TM mode is adjusted by the horizontal slot in this structure, realizing an ultra-flat dispersion profile, which varies from −10.6 ps·nm−1·km−1 to +11.14 ps·nm−1·km−1 over a wavelength range of 1080 nm (from 1760 nm to 2840 nm).

2. Waveguide design and dispersion tailoring by the slot structure

Figure 1(a) shows the reversed-rib waveguide structure. Due to the low transformation temperature of As2S3,[28,31] the fabrication process of this structure is relatively simple. Firstly, a groove was made on the silica substrate by wet etching or dry etching. Secondly, the As2S3 film was deposited by traditional vacuum deposition methods (thermal evaporation, chemical vapor deposition or sputtering). Then, the sample was annealed to melt the As2S3 material, which made the upper surface of the film flat to realize the reversed-rib structure. It can also be fabricated by the micro-trench filling method based on solution-processed ChG.[27,29,51] Finally, a silica layer as the cladding of the waveguide could be fabricated by sputtering or the method based on the silicon based polymer. In this process, there is no photolithography or etch processes on the As2S3 film, which reduces the difficulty in the fabrication of ChG waveguides. It is worth noting that this structure can also be fabricated by the traditional top-down method, in which the ChG glass and silica materials are deposited layer by layer firstly, and then only the last ChG layer needs to be etched to form the ridge structure.

Fig. 1. Sketches of the proposed waveguide structures. (a) The reversed-rib waveguide without slot. (b) The reversed-rib waveguide with a silica slot.

To tailor the dispersion of the reversed-rib waveguide, we introduce a horizon slot above the rib, which is shown in Fig. 1(b). It can be realized by simply repeating the deposition of both As2S3 film and silica layer one more time.

Dispersion tailoring is carried out by properly optimizing the structural parameters of the waveguide including the rib width (W1), the rib height (H1), and the thicknesses of the As2S3 layer on the subtract (H2), the silica slot (Hs), and the As2S3 layer above the slot (H3). In order to prevent the oxidation of the ChG glass, the waveguide is covered with a layer of SiO2 with a thickness of Hcladding. In this paper, the effective indexes (neff) of the fundamental quasi-TE mode and quasi-TM mode with different wavelengths are calculated by the finite-element mode method. In the calculation, the material dispersions of As2S3[52] and silica[53] are taken into account using their Sellmeier equations. The dispersion can be calculated by

where c and λ are the speed of light and the wavelength in vacuum, respectively. Figure 2 shows the dispersions of the fundamental quasi-TE/TM modes in the waveguides shown in Fig. 1. The dashed-dotted blue line in Fig. 2 is the material dispersion of the As2S3. The red and black dashed lines are the dispersion of the fundamental quasi-TE and TM modes in the reversed-rib waveguide without slot (as shown in Fig. 1(a)). Here, the width and height of the rib are W1 = 1400 nm and H1 = 714 nm, respectively. The thicknesses of the As2S3 layer above the substrate (H2) and the silica cladding layer (Hcladding) are 80 nm and 180 nm, respectively. Calculations show that no higher-order mode is supported under these structure parameters, except for the fundamental quasi-TE/TM modes. It can be seen that both modes have parabola-like dispersion curves with two zero dispersion wavelengths (ZDWs). The maximums of the dispersions of the quasi-TE and TM mode are close to 200 ps·nm−1·km−1 and 320 ps·nm−1·km−1, respectively. Besides, it can be seen that the near-zero dispersion bands are all narrow near the ZDWs of the two modes in this waveguide.

Fig. 2. Dispersion of the fundamental quasi-TE and quasi-TM modes in the reversed-rib waveguide with/without the slot. The solid red and black lines are the dispersion of the quasi-TE and quasi-TM modes in the reversed-rib waveguide with a slot, respectively. The dashed red and black lines are the dispersion of the quasi-TE and quasi-TM modes in the reversed-rib waveguide without slot, respectively. The blue dashed-dotted line is the material dispersion of the As2S3.

The red and black solid lines in Fig. 2 are the dispersion of the fundamental quasi-TE and -TM modes in the reversed-rib waveguide with the slot (as shown in Fig. 1(b)). Here, the width and height of the rib (W1, H1), the thickness of the As2S3 layer above the substrate (H2), and the silica cladding layer (Hcladding) are all the same as those of the waveguide without the slot. The thicknesses of the silica slot (Hs) and the As2S3 layer above the slot (H3) are 93 nm and 230 nm, respectively. It can be seen that, for both modes, the dispersion shows flatter profiles with broader near-zero dispersion bands than that of the reversed-rib waveguide without the slot. The inset of Fig. 2 shows the detail of the dispersion of the quasi-TM mode in the reversed-rib waveguide with the slot. It can be seen that it has two ZDWs at 1849 nm and 2742 nm, respectively. The dispersion is limited between −10.6 ps·nm−1·km−1 and +11.14 ps·nm−1·km−1 over a band of 1080 nm (from 1760 nm to 2840 nm).

To show the mechanism of the dispersion tailoring effect of the silica slot, the electric field distributions of the quasi-TM modes with different wavelengths in the proposed waveguide with slots and without slots are calculated and shown in Figs. 3(a) and 3(b), respectively. It can be seen that the electric field concentrates in the rib region at the short wavelength. As the wavelength increases, the electric field gradually transfers into the SiO2 slot. This mode field transfer provides a modification on the waveguide dispersion of the reversed-rib waveguide. Hence, the slot provides additional dispersion tailoring freedom.

Fig. 3. Quasi-TM mode evolution of the waveguide for various wavelengths (a) with slots and (b) without slots.
3. The influence of structure parameters on the dispersion property

The slot structure provides enough degrees of freedom to tailor the waveguide dispersion. On the other hand, the variation of structure parameters, which is unavoidable in waveguide fabrication, would also impact the dispersion property of the waveguide. Figure 4 shows the influence of structural parameters on the dispersion of the fundamental quasi-TM mode by changing them around the optimized values.

Fig. 4. Dispersion profiles of the slot waveguides with different (a) slot thicknesses (Hs), (b) heights of groove (H1), (c) heights of the As2S3 layer (H2), (d) widths of groove (W1), and (e) heights of another As2S3 layer (H3).

Figure 4(a) shows the calculation results of dispersion with different thicknesses of the slot. The other structure parameters used in this calculation are the same as the optimized design used in Fig. 2 (This condition is also used in the calculations of Figs. 4(b)4(e)). It can be seen that the dispersion curve with a larger slot thickness shows higher suppression at the long-wavelength band. By optimizing the slot parameters, a dispersion curve with a flat profile can be realized.

The dispersion profiles for different rib heights (H1) are plotted in Fig. 4(b). It can be seen that the dispersion profile is almost unchanged. While it rises as a whole as H1 increases from 704 nm to 724 nm, the maximum dispersion value rises from 0 ps·nm−1·km−1 to 21.43 ps·nm−1·km−1, with a rate about 1 ps·nm−1·km−1 per nanometer. A similar effect can be seen when the thickness of the As2S3 layer on the subtract (H2) changes. Figure 4(d) shows the calculation results under different rib widths (W1). When the width varies from 1390 nm to 1410 nm, the dispersion shape and slope at the short wavelength are almost unchanged. However, the dispersion in the longer wavelength region changes rapidly and the dispersion on the peak increases from 9.5 ps·nm−1·km−1 to 12.4 ps·nm−1·km−1 with a rate about 0.15 ps·nm−1·km−1 per nanometer.

The influence of the thickness of the As2S3 layer above the slot (H3) is calculated and shown in Fig. 4(e). The dispersion profile drops with the increase of H3. The maximum dispersion reduces with a rate of about 1 ps·nm−1·km−1 per nanometer and shifts to longer wavelength.

By comparing the results in Fig. 4, it can be seen that the slot thickness is the most sensitive structure parameter for the dispersion profile. Hence, its fine control in the fabrication process is the key to obtain a flat and near-zero dispersion profile.

4. The performance of the proposed waveguide as a third-order nonlinear waveguide

The As2S3 glass has the third-order nonlinear index (n2) of 2.92 × 10−18 m2/W and a low TPA coefficient of 6.2 × 10−15 m/W at the telecom band. Its nonlinear FOM is as high as 312, which is much higher than that of the silicon waveguide (0.77).[31] On the other hand, the proposed waveguide structure provides an effective way to realize a flat dispersion profile, which is preferred in many third-order parametric nonlinear processes. Hence, the proposed waveguide has great potential as a nonlinear medium for developing integrated nonlinear devices. To show this more clearly, we calculated the third-order nonlinear coefficient of the fundamental quasi-TM modes in the proposed waveguide with different wavelengths by

where F(x,y) is the distribution of the electric field and is the light wavelength, Aeff is the effective area of the mode, which is calculated by[54]

The structure parameters used here are the same as those used in Fig. 2, in which quasi-TM mode of the proposed waveguide has an ultra-flat profile with a broad near zero dispersion band. Figure 5 shows the calculated (blue circles) and (black squares), respectively. It can be seen that, as the wavelength increases, the effective area of the mode increases while the nonlinear coefficient reduces. At a short wavelength, most of the electric field is confined in the reversed-rib of the As2S3 waveguide. As the wavelength increases, the electric field increasingly extends to the silica slot. Hence, the nonlinear coefficient gradually decreases due to the electric field extension into the silica slot with a lower nonlinear index.

Fig. 5. The effective area Aeff (blue circles) and the nonlinear coefficient γe (black squares) of the fundamental quasi-TM mode in the proposed reversed-rib waveguide with a slot. The structure parameters used in the calculation are the same as those used in Fig. 2.

To demonstrate the potential of this waveguide in the third-order parametric processes, we calculated the linear phase mismatching of the degenerate four-wave mixing by Δβ = βs + βi − 2βp, where βs, βi, and βp are propagation constants of the signal, idler, and pump lights, respectively. The quasi-TM mode with the ultra-flat profile is considered in this calculation. The calculation results are shown in Fig. 6. Figure 6(a) is the linear phase mismatching under different pump wavelengths λpump and signal-pump wavelength differences Δλ. Different colors in the figure represent the values of linear phase mismatching under different λpump and differences Δλ. It can be seen that if the pump light is set close to the two ZDWs, small linear phase mismatching could be realized over large regions of signal-pump wavelength difference Δλ.

Fig. 6. Calculation results of linear phase mismatching of the degenerate FWM process in this waveguide. (a) Linear phase mismatching Δβ under different pump wavelengths and signal-pump wavelength differences Δλ, (b) linear phase mismatching under different signal wavelengths when the pump wavelength is set near the shorter ZDW, (c) linear phase mismatching under different signal wavelengths when the pump wavelength is set near the longer ZDW.

To show this more clearly, the linear phase mismatch profiles when the pump wavelength is set near the shorter ZDW (1849 nm) and longer ZDW (2742 nm) are plotted in Figs. 6(b) and 6(c), respectively. It can be seen that when pump wavelength is set to 1853 nm, which is close to the shorter ZDW (1849 nm), the linear phase mismatching profile is like an upward parabola with a zero point at Δλ = 0 and a broad near-zero band. When the pump wavelength decreases, e.g., to 1833 nm, the parabola-like profile shrinks, leading to a smaller near-zero band. On the other hand, when the pump wavelength increases, e.g., to 1863 nm, two regions with negative linear phase mismatching appear at two sides, which would be deeper and wider if the pump wavelength increases further. The profile of the linear phase mismatching when the pump light is near the longer ZDW (2742 nm) is similar with the case that it is near the shorter ZDW. However, its variation with the pump wavelength is different. It can be seen that the two regions with negative linear phase mismatching appear when the pump wavelength is shorter than the longer ZDW, with a broader near-zero band.

It is well known that a small negative linear mismatching is preferred to the degenerate fourwave mixing because it can be compensated by the nonlinear mismatching term under a proper pump level. Hence, to realize the broad band degenerate fourwave mixing in the quasi-TM mode of the proposed waveguide, the pump light should have a wavelength a little longer than the shorter ZDW or a little shorter than the longer ZDW.

5. Conclusion

In this paper, we proposed a reversed-rib As2S3 waveguide with a slot to tailor its dispersion. By properly designing its structure parameters, an ultra-flat and near zero dispersion profile can be realized, in which the dispersion is limited between −10.6 ps·nm−1·km−1 and +11.14 ps·nm−1·km−1 over a band of 1080 nm. The slot structure provides a dispersion tailing effect to extend the low dispersion band, which is due to the electric field transfer process from the reversed-rib waveguide to the slot as the wavelength increases. To demonstrate its potential as a nonlinear waveguide for integrated nonlinear optical devices, its nonlinear coefficient and the linear phase mismatching profile for the degenerate FWM are also calculated, showing that the waveguide can support broadband degenerate FWM in near infrared and middle infrared band if the pump wavelength is set a little longer than the shorter ZDW or a little shorter than the longer ZDW. The proposed As2S3 waveguide can be fabricated without photolithography and etching processes on the As2S3 film, which can reduce the difficulty in the fabrication of ChG waveguides, showing great potential in applications of nonlinear photonic devices.

Reference
1Chavez B J MBodenmüller DFremberg THaynes RRoth M MEisermann RLisker MZimmermann LBöhm M 2014 J. Opt. Soc. Am. 31 2846
2Lamont M RLuther-Davies BChoi D YMadden SEggleton B J 2008 Opt. Express 16 14938
3De Leonardis FPassaro V M N 2011 Adv. Optoelectron. 2011 751498
4Duan LYang Z YLiu CYang W L 2016 Chin. Phys. Lett. 33 10501
5Wang H YXu W CLuo Z CLuo A PCao W JDong J LWang L Y 2011 Chin. Phys. Lett. 28 024207
6Amiri I SAfroozeh ABahadoran MAmiri I SAfroozeh ABahadoran M 2011 Chin. Phys. Lett. 28 104205
7McCarthy JBookey HBeecher SLamb RElder IKar A K 2013 Appl. Phys. Lett. 103 151103
8Zhang HDas SHuang YLi CChen SZhou HYu MGuo-Qiang L PThong J T L 2012 Appl. Phys. Lett. 101 021105
9Yu YGai XWang TMa PWang RYang ZChoi D YMadden SLuther-Davies B 2013 Opt. Mater. Express 3 1075
10Zhang LLin QYue YYan YBeausoleil R GWillner A E 2012 Opt. Express 20 1685
11Liang T KTsang H K 2004 Appl. Phys. Lett. 84 2745
12Liu QGao SLi ZXie YHe S 2011 Appl. Opt. 50 1260
13Lin QZhang JFauchet P MAgrawal G P 2006 Opt. Express 14 4786
14Lamont M Rde Sterke C MEggleton B J 2007 Opt. Express 15 9458
15An LLiu HSun QHuang NWang Z 2014 Appl. Opt. 53 4886
16Jin BYuan JYu CSang XWei SZhang XWu QFarrell G 2014 Opt. Express 22 6257
17Collins M JClark A SHe JChoi D YWilliams R JJudge A CMadden S JWithford M JSteel M JLuther-Davies BXiong CEggleton B J 2012 Opt. Lett. 37 3393
18Cardinal TRichardson K AShim HSchulte ABeatty RLe Foulgoc KMeneghini CViens J FVilleneuve A 1999 J. Non-Cryst. Solids 256 353
19Viens J FMeneghini CVilleneuve AGalstian T VKnystautas E JDuguay M ARichardson K ACardinal T 1999 J. Light. Technol. 17 1184
20Feigel AKotler ZSfez BArsh AKlebanov MLyubin V 2000 Appl. Phys. Lett. 77 3221
21Zhai YQi RYuan CDong SZhang WHuang Y 2016 IEEE Photonics 8 2700709
22Al-Kadry ALi LAmraoui M ENorth TMessaddeq YRochette M 2015 Opt. Lett. 40 4687
23Zhai YYuan CQi RZhang WHuang Y 2015 IEEE Photonics 7 7801609
24Chiles JMalinowski MRao ANovak SRichardson KFathpour S 2015 Appl. Phys. Lett. 106 111110
25Luther-Davies B Yu Y Zhang B Gai X Zhai C Qi S Guo W Yang Z Wang R Choi D Y Madden S Moller U Kubat I Petersen C Brilland L M’echin D Caillaud C Troles J Bang O 2015 Nonlinear Optics 2015, OSA Technical Digestp. NTu1A.410.1364/NLO.2015.NTu1A.4
26Yu YGai XMa PChoi D YYang ZWang RDebbarma SMadden S JLuther-Davies B 2014 Laser Photonics Rev. 8 792
27Zou YMoreel LLin HZhou JLi LDanto SMusgraves J DKoontz ERichardson KDobson K DBirkmire RHu J 2014 Adv. Opt. Mater. 2 759
28Zou YZhang DLin HLi LMoreel LZhou JDu QOgbuu ODanto SMusgraves J DRichardson KDobson K DBirkmire RHu J 2014 Adv. Opt. Mater. 2 478
29Zha YLin P TKimerling LAgarwal AArnold C B 2014 ACS Photonics 1 153
30Gai XMadden SChoi D YBulla DLuther-Davies B 2010 Opt. Express 18 18866
31Eggleton B JLuther-Davies BRichardson K 2011 Nat. Photonics 5 141
32Eggleton B J 2010 Opt. Express 18 26632
33Karim M RRahman B M AAzabi Y OAgrawal AAgrawal G P 2015 JOSA B 32 2343
34Luther-Davies BGai XMadden S JChoi D YYang ZWang RMa P Yu I 2013 CLEO: Science and Innovations (Optical Society of America)p. CM1L
35Eggleton B JVo T DPant RSchr JPelusi M DYong C DMadden S JLuther-Davies B 2012 Laser Photonics Rev. 6 97
36Collins M JClark AHe JShahnia SWilliams R JJudge A CMagi EChoi D YLuther-Davies B Eggleton B J 2012 Frontiers in Optics (Optical Society of America)p. FTu4D
37Qiao H AAnheier N CMusgrave J DRichardson K Hewak D W 2011 SPIE Defense Security, and Sensing (International Society for Optics and Photonics)p. 80160F
38Bao C Yan Y Zhang L Yue Y Ahmed N Agarwal A M Kimerling L C Michel J Willner A E 2015 J. Opt. Soc. Am. 32 26
39Liu YYan JHan G 2014 Appl. Opt. 53 6302
40Wang SHu JGuo HZeng X 2013 Opt. Express 21 3067
41Nolte P WBohley CSchilling J 2013 Opt. Express 21 1741
42Zhu MLiu HLi XHuang NSun QWen JWang Z 2012 Opt. Express 20 15899
43Zhang LYue YXiao-Li YWang JBeausoleil R GWillner A E 2010 Opt. Express 18 13187
44He JXiong CClark A SCollins M JGai XChoi D YMadden S JLuther-Davies BEggleton B J 2012 J. Appl. Phys. 112 123101
45Choi D YMaden SRode AWang RLuther-Davies B 2008 J. Non-Cryst. Solids 354 3179
46Gai XHan TPrasad AMadden SChoi D YWang RBulla DLuther-Davies B 2010 Opt. Express 18 26635
47Choi D YMadden SRode AWang RBulla DLuther-Davies B 2008 J. Non-Cryst. Solids 354 5253
48Choi D YMadden SBulla D AWang RRode ALuther-Davies B 2010 IEEE Photonics Technol. Lett. 22 495
49Pittman T BJacobs B CFranson J D 2005 Opt. Commun. 246 545
50Zhai YQi RYuan CZhang WHuang Y 2016 Appl. Phys. Express 9 052201
51Barile C JNuzzo R GGewirth A A 2015 J. Phys. Chem. 119 13524
52Rodney W Set al.1958Anon Refractive Index of As2S3 (Arsenic Trisulfide)
53Malitson I H1965Anon Refractive Index of SiO2 (Silicon Dioxide, Silica, Quartz)
54Li LZou YLin HHu JSun XFeng N NDanto SRichardson KGu THaney M 2013 J. Light. Technol. 31 4080